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Abstract 

The rotational properties of an attractively interacting Bose gas are studied using analytical and 
numerical methods. We study perturbatively the ground state phase space for weak interactions, 
and find that in an anharmonic trap the rotational ground states are vortex or center of mass 
rotational states; the crossover line separating these two phases is calculated. We further show 
that the Gross-Pitaevskii equation is a valid description of such a gas in the rotating frame and 
calculate numerically the phase space structure using this equation. It is found that the transition 
between vortex and center of mass rotation is gradual; furthermore the perturbative approach is 
valid only in an exceedingly small portion of phase space. We also present an intuitive picture of 
the physics involved in terms of correlated successive measurements for the center of mass state. 
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I. INTRODUCTION 



The rotational properties of atomic Bose condensed gases have been a matter of intense 
scientific interest recently. The main role in these studies is played by vortices, which are the 
manifestations of quantized circulation in the quantum gas in question. However, vortices 
are by no means the only possible rotational states. 

By rotating a condensate in a harmonic trap with a certain frequency Q a vortex can 
be created; when rotating the condensate a bit faster a second vortex appears, and so 
on . Lattices with a large number of vortices have been experimentally observed in many 
groups . There is an upper limit, however, for the rotation frequency in a harmonic 
trap. This limit is the trap frequency u of the confining trap, above which the center of 
mass of the condensate is destabilized ^ . Interesting physics has been done by the Boulder 
group working only slightly below this limit P]. Having prepared a condensate with a vortex 
lattice, atoms with low angular momenta are then removed. Reaching the rotation rate 0.95u; 
a metastable giant vortex with a phase singularity on the order of 50 is observed. 

On the other hand, the trap limit can be crossed if the trap is not purely harmonic but 
the potential has a quartic term also j^. The steeper potential makes fast rotation possible 
and thus gives a tool to test theoretical predictions including multiply quantized vortices ^ , 
strongly correlated vortex liquids or Quantum Hall like states j?,!^. The first experimental 
studies have already been done at ENS 

The rotation of a condensate in an anharmonic trap is especially interesting when the 
effective interactions of the atoms are attractive. It was shown by Wilkin et al. that for 
weak attractive interactions and a given angular momentum, in the lowest energy state the 
motion is carried by the center of mass (c. m.) [3]. For this state the single-particle reduced 
density matrix has more than one macroscopic eigenvalue in the laboratory reference of 
frame. This might indicate that the cm. rotational state is a fragmented condensate, but 

n 

actually, the state is Bose condensed if it is viewed in the c. m. frame Unfortunately, in 
a harmonic-oscillator trap this state is never thermodynamically stable because the critical 
frequency for the excitation is equal to the frequency where the gas is destabilized. An 
anharmonic term in the trapping potential removes the instability and makes the state of 
c. m. rotation, as well as vortex states, attainable 12]. 

In this paper we shall study the rotational states of an attractive zero temperature Bose 
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gas confined in an anharmonic trap. We begin by discussing the nature of tlie relevant states 
of motion and perturbation tlieory in Sec. |n] In Sec. II I II we discuss the validity of different 
equations of motion. The problem is approached within the ordinary Gross-Pitaevskii (GP) 
mean field theory in Sec. IIVI Section |3 deals with the correlation properties of the c. m. 
rotational state at the limit of small anharmonicity. Concluding remarks are made in Sec. IVll 



II. HAMILTONIAN AND PERTURBATION THEORY 



The Hamiltonian for atoms with binary s-wave interactions is 

t2 AT .2 N jj 

^--I^Ell + EVW + ^E^lr.-r,) (1) 

1=1 « 1=1 l^J 

where the strength of the contact interaction Uq = Anh^a/m with the scattering length a 
which is here assumed negative. The external field V consists of a cylindrically symmetric 
harmonic potential with a small quartic addition in the radial direction: 



V{r, 6, z) = -m 
2 



"osc 
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(2) 



Here the dimensionless parameter A describes the strength of the anharmonic term and the 
oscillator length a^sc is defined as Oosc = {h/mcuY^'^. 

In Refs. [3, 1^, the harmonically trapped gas (A = 0) was studied by treating the interac- 
tion energy perturbatively, in which case the many-body eigenstates are products of nodeless 
harmonic-oscillator single-particle eigenstates. Here we shall generalize this approach by in- 
cluding also the anharmonic potential 2Q\. To this end, we first pass to dimensionless units, 
where lengths are scaled by Oosc and energies by hu. We obtain 

^ = E i + ^ E + + ^ E ^(r. - r,), (3) 

i=l i 1=1 j^j 

where we have defined an effective 2D interaction g by integrating out the 2;-component 



wave function to obtain g = y uJz/{27Tuj)a/aosc (this can be done because the radial-axial 
decoupling is exact in perturbation theory). In the attractive case, g is negative. We shall 
retain the symbols r, H etc. for the dimensionless quantities in this section, because no 
confusion can arise. 

A general form for a many-body state in two dimensions is 

N 

(4) 

mi,...,mjv j=l 



where 0mn(r) is the harmonic-oscillator single-particle eigenstate in two dimensions with 
angular momentum m and radial quantum number n; only states with n = participate 
in perturbation theory. In the unperturbed case (A = a = 0), for a given total angular 
momentum L all states which fulfill mi + . . . + m^r = L with all mj > (and have no 
radial nodes, i. e. n = 0) are degenerate. Thus we have, in principle, to perform degenerate 
perturbation theory on a vast subspace of Hilbert space (the number of basis states is 
(A^ + L — 1)!/[A^!(L — 1)!]). However, from physical arguments and known exact results, we 
obtain guidance to what the relevant candidate states should be. 

In the purely harmonic case with finite attractive interactions, the ground-state mam 
body wave function with angular momentum L and particle number N was found to be [iC 



= expJ-f:|.,|V2j, (5) 

where zj = xj + iyj are the particle coordinates. The quantity within parentheses is the 
c. m. coordinate; therefore, this wave function describes rotation of the c. m. as already 
mentioned in the introductory paragraph. We consider this state as the ideal c. m. rotational 
state. 

When instead the anharmonicity A is nonzero and the interactions vanish, it was shown 
in Ref. 0] that the ground state is the vortex state 

N 1 

^l{zu...,zr,) = l[^^^T^-^''^'^^ (6) 

which is just a Bose-Einstein condensate containing an m-fold quantized vortex, thus having 
angular momentum Nm. Vortex-array states do not come into question because they have 
a lower density on the average, and are therefore not favored by the attractive interaction. 

It should be noted that Eqs. (0) and (jUj) are the extreme cases where either the interaction 
or anharmonicity are absent. In a more general context we make the difference between these 
states in terms of cylindrical symmetry. The vortex state is cylindrically symmetric and the 
ideal c. m. rotational state corresponds clearly to the nonsymmetric case. Between these 
two states, however, there is a region of nonsymmetric cases, where the c. m. is in rotation 
although these states do not correspond to the ideal c. m. rotational case of Eq. (0). The 
eading instability of the vortex state towards a nonsymmetric state was examined in Ref. 
We include these cylindrically nonsymmetric states into our general definition of c. m. 



rotational states. Other interpretations are possible, though, see 



li 



There now remains to compare the energies of the two states (jSJ and © 21|. The 
interaction energies of the two states were calculated in Ref. and amount to 



{^1\H,^,\^1) = gN{N-l 

) 

(2m- 1)!! 



The absolute magnitude of the interaction energy in the vortex state is smaller than that 
of the ideal c. m. state, and therefore the latter is favored by the attractive interactions 
(remember that g <Q). The quartic energies, on the other hand, are proportional to 

(^flVr^l^f) = ivf2 + 4- + — ---V 

(v^:;i5:rj|vl/:;) = iV(m + l)(m + 2), (8) 

so that for equal angular momentum, the quartic energy of the vortex state is the lowest. 

In the vortex state the angular momentum per particle m can only take on integer values. 
In contrast, in the ideal c. m. state, q = L/N is only quantized in quanta of so when 
we pass to the infinite limit g is a continuous variable. At fixed L, it is easy to determine 
the boundary between the rotational phases: the ideal c. m. state is energetically favorable 
when 

A (2m-l)!!\ 

-m < ^ 1 - ^ — ^ . 9) 

We have discarded all terms 0{N~^). For m = q = 1 the condition for c. m. rotation 
is A < \g\N and for faster rotation, the regime where the ideal c. m. state is favorable is 
smaller. Since the right-hand side of Eq. © increases slower than the left-hand side as m 
increases, we conclude that for any given pair of parameters {X,gN), there exists a smallest 
angular momentum L above which a vortex state is favorable. 

We now change variables and work at fixed angular velocity Q rather than fixed L. In 
the ideal c. m. state we can simply differentiate the energy: 

n = ^ = 1 + A(2 + g) ^ g = - 2. (10) 

uL/ A 

The critical frequency for the excitation of rotational motion is thus = 1 + 2A, and 
above this threshold we substitute the above expression for q in order to obtain the energy 
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function of Vl: 

1 (Q-lV 

-[Ec{n, A, g) - QL] = -l + 2fi - ^ ^ - A - |^7|iV. (11) 

The energy of the vortex state as a function of m is 

1 1 (2m — IV 

-E,(m) = 1 + m + -A(m + l)(m + 2) - \g\N^—^—^. (12) 

Comparing Eqs. (fTTj) and (fT^. the critical angular velocities for the successive vortex states 
are found to be 

n^ci = l + 2\ + ^\g\N, 
= l + 3X + ^\g\N, 

(2m — 3V 

^l,m = 1 + (m + 1)A + m^-^^-^. (13) 

For given gN, A and Q it can now be determined which of the states and has the 
lower energy. When < Q < fi^2' ^^"^ that the ideal c. m. state is favorable if the 
following equivalent conditions hold: 

(Q - 1)^ 

\g\N > 6{n - 1) - ^ — - 8A, 

A 



A < |(n - 1) - IwA' - Jl(a - ir - l-^m(n - 1) + ^^m, 



n < l + 3A-A/A(A-|^|iV). (14) 



If A- y/A(A - \g\N) < \g\N/2, then the vortex state is favorable for all > i^^^. If \g\N > A, 
then the ideal c. m. state is the favorable one in the whole interval Q'^i < Q < but 
the vortex state may be favorable for higher frequencies, as discussed above. To summarize 
the dependencies: (1) a small anharmonic term A favors the ideal c. m.-rotation state; (2) 
a small interaction term \g\N favors the vortex state; (3) a large angular velocity Q favors 
the vortex state. 



III. EQUATIONS OF MOTION 

The perturbative approach has, naturally, a limited range of validity; as we shall see, it 
is in fact accurate only for very small values of gN and A. There is thus need for a more 
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Vnir) + t/o|<f (r)p| ^(r) = ^$(r). (16) 



general scheme and we shall now discuss what kind of approximation can be used to describe 
the attractive gas. 

In a purely harmonic potential the center of mass motion decouples from the internal 
motion. This is not the case in an anharmonic trap, but there is still an approximate 
decoupling for weak anharmonicity. By taking advantage of this, and the fact that the 
internal motion is Bose condensed, coupled equations of motion for the c. m. and internal 
motion were derived in Ref. jl^ : 

{-^^2 + ^c(R)| ^c(R) = EcM^), (15) 
\ 2m dv"^ 

Here, the c. m. is described quantum mechanically by the wave function ipci^)-, and the 
Bose- Einstein condensed internal motion is governed by the condensate wave function $(r), 
where R is the c. m. coordinate and r is the particle coordinate relative to the c. m. The 
effective potentials for the center-of-mass and relative motion are 

Vc{R) = Imu;\1 + 4X^^)R' + (17) 

^ ^osc ^'^osc 

Vn{r) = -mu\l + 4\^y + -y^r\ (18) 

^ '^osc ^'^osc 

where (i?^) = / A^^c?Ri?^|^/'c(R)P and (r^) = J (irr^|$(r)p gives the coupling between 
the two wave functions. These equations are capable of describing both the (ideal) state of 
c. m. rotation and vortex states and were used in Ref. to find the parameter regimes for 
those two states of motion. The perturbative results of Sec. [H] emerge from these equations 
in the limit A — > 0, gN — > if is assumed large. 

Here, we shall take an alternative path to describing the motion of the attractive gases. 
Namely, we argue that the ordinary Gross-Pitaevskii equation_ 15[ can be used to describe 



these rotating systems. The reason is simple: as shown in Ref. pj|, the state of c. m. rotation 
is Bose-Einstein condensed according to an observer co-moving with the c. m. Therefore, 
if we transform to a coordinate system moving with the c. m., we should be able to use 
the Gross-Pitaevskii equation. But if the cloud is in its ground state, under a rotational 
force with frequency Q, the c. m. is also rotating with the frequency Q and the coordinate 
transformation is effected simply by going to a rotating frame as usual : 

H ^ H' = H -VlL, (19) 



where L is the angular momentum operator. Solving the Gross-Pitaevskii equation in the 
rotating frame should thus allow us to describe stationary states of the cloud, both in the 
vortex and the c. m. states. This result comes as a bit of a surprise, considering that in 
Ref. Jicl] the c. m. state was found not to be Bose-Einstein condensed. We note also that 
working at fixed Vt is the correct description of the common experimental approach of letting 
the cloud equilibrate under the influence of a rotational drive. 

Compared with the two coupled wave equations (|15II18|) . the full Gross-Pitaevskii equation 
is more general: it lifts the restriction of small values of A, so that rotational motion in any 
trap can be described as long as one can move to a reference frame in which the motion is 
Bose-Einstein condensed. On the other hand, the GPE approach describes the c. m. motion 
classically instead of quantum mechanically, but this is hardly an issue unless one wishes to 
specifically study the quantum fluctuations in the c. m. motion. The quantum mechanical 
c. m. motion will be illustrated in Sec. 



IV. NUMERICAL SIMULATIONS 



As we argued in the preceeding section, the attractive Bose gas is well described by the 
single component Gross-Pitaevskii equation in a frame rotating with frequency f2. 



___V2 + V{v)-Vt-L + Uo\^ (r) 



^ (r) = /i^ (r) . (20) 



To solve this numerically for the anharmonic potential of Eq. Q, we choose A = 0.15, 
= 1000, = 27r X 30 Hz, a;^ = 27r x 180 Hz and the mass of atomic Li^. For the harmonic 
trap frequency much stronger in z-direction, we may assume the motion in this direction to 
be frozen in the ground state of the trap and the problem becomes two-dimensional. The 
ground state is now found by numerically propagating the time-dependent counterpart of 
Eq. ()20|) in imaginary time. The grid size varies from 128 x 128 points up to 256 x 256. 

With the fixed parameters mentioned above, we map the ground state phase-space with 
respect to the rotational frequency and the scattering length. The resulting phase-space 
diagram is plotted in Fig. [T] The result of the decoupling approximation, Eqs. ()15lll8j) . 



with a Gaussian ansatz for the density 



12j | is also included. We see that the ground state 



is nonrotating below Vt/uo ^ 1.2. The cloud just stays at the bottom of the trap with 
angular momentum L = (within the numerical uncertainty). As the rotation frequency is 



8 



increased the ground state depends on the interaction strength. For weak interactions the 
angular momentum L/N is quantized in integer values that increase with rotation frequency. 
The ground state is a (multiply) quantized vortex (see Fig. |2)). On the contrary, stronger 
interactions break the rotational symmetry. The circulation is not quantized anymore and 
the density of the condensate has a shape of a crescent as can be seen from Fig. |21 This 
configuration corresponds to the c. m. rotational state. For even stronger interactions, the 
density becomes more concentrated and the cloud attains an approximately ellipsoidal shape 
residing off center. 

Interestingly, the coupled equations of motion ()15II18|1 . derived from the factorization of 
the wave function into c. m. and internal parts, fail to describe the elongated shapes that 
the cloud attains for moderately strong attraction. The assumptions behind those equations 
imply, namely, cylindrical symmetry of the condensate wave function. The factorization 
ansatz thus fails precisely in the situation when the quartic potential has a strong effect on 
the shape of the cloud, i.e., when we are close to the phase boundary. 

From the perturbative method we have an expression for q = L/N in Eq. as a 
function of Q and A. We perform a comparison between this simple perturbative formula 
for the ideal c. m. state and the GP procedure by choosing the set of these pairs {Q, A) to 
be such that q = 1- In Fig. El we plot the results from the GP simulations for the angular 
momentum per atom as a function of A. When the anharmonicity gets smaller q seems to 
approach unity. Still, for the given parameter range the difference between numerical and 
perturbative results is clear. This is somewhat unexpected, so to verify the numerical results 
we also plot the corresponding values of q for Gaussian trial wave functions in the coupled 
Eqs. ()15II18|) . The results are comparable to the ones from GP simulations; this shows that 
perturbation theory can only be trusted in an extremely small portion of phase space. 

That the anharmonic trapping potential is essential for stabilizing these states of motion 
can be seen from the following argument. The effective potential seen by the rotating atoms 
is the sum of the actual potential V{r, 9, z) and the centrifugal term: 

KfT(r) = \m{uj^ - n'^y + \mu:^\^. (21) 

For slow rotation, Vt < uj, the prefactor of the term is positive and the effective potential 
is just a more shallow anharmonic potential. In that regime there is no rotation for the 
attractive gas. For rotation faster than the trap frequency, however, the effective potential 
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FIG. 1: The ground state phase-space diagram of a condensate in a rotating trap with A = 0.15 as 
a function of the rotating frequency and the scattering length. Different ground states are marked 
by the following symbols: Q, nonrotating state; V) vortex state (the number is the circulation of 
a multiply quantized vortex); *, c. m. rotational state. The solid lines give the same phase-space 
diagram obtained from the equations of motion where the c. m. and internal motion are decoupled, 
using a Gaussian ansatz for the wave functions. 



FIG. 2: Density plots for different ground state configurations. Bright shades indicate high density. 
On the left, we have a doubly quantized vortex for a = —1.0 a.u., in the center c. m. rotational 

state for a = —4.0 a.u. and on the right c. m. state for a = —10.0 a.u. Anharmonicity is A = 0.15 

1 

and the rotation frequency is fixed to Q/u = 1.35. The unit of length for x and y is {h/mu) 2 . 
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FIG. 3: The angular momentum (in units of K) per atom as a function of the dimensionless anhar- 
monicity parameter A for the c. m. state. The parameter JI/lj is chosen such that in perturbation 
theory LjN = 1. Because we compare this to the results from GP simulations (marked by the sym- 
bol *) we have to fix the dimensionless interaction parameter gN also. The parameter set for the 
lowest A is (A = 0.01; Vi = 1.03; gN = 0.015) and the others are (0.02; 1.06; 0.03) , (0.05; 1.15; 0.075) 
and (0.1; 1.3; 0.15) respectively. The dashed line is variationally calculated for Gaussian trial wave 
functions in the decoupling approximation. 

attains a Mexican-hat shape and the ground-state density distribution lies along the bottom 
of this toroidal potential. If the trap bottom is deep and the effective trap width is small 
the condensate is confined in an effectively one- dimensional torus which has been studied 
analytically in two recent articles 3,^^- Just as in the present study, two types of state 
are found in this idealized geometry, termed the uniform- density state and the localized 
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FIG. 4: Density plot of a c. m. rotational state obtained by solving 3D GP equation 
z = plane and on the right y = plane. The relevant parameters are 0,/uj = 1.40, a 
and A = 0.15. The unit of length for x, y and z is (h/muj)^ . 



. On the left 
= -10.0 a.u. 



density state (bright soliton). Clearly, the uniform- density state can be identified with our 
vortex state, and the bright soliton with the c. m. state. The conclusion drawn from the 
GP analysis is thus that the crossover from vortex to cm. rotation is gradual, whereas the 
coupled equations (jl5lll6|) predicted a discontinuous crossover |12]. In view of the arguments 
presented in Sec. IIIII to the effect that the GP equation provides a valid description of the 
gas, the correct conclusion is that the crossover is gradual. This conclusion is supported by 
the exact findings in the simplified ID model Q, 17 1. The failure of the coupled-equation 
approach can be ascribed to the fact that it is restricted to cylindrically symmetric density 
profiles and is thus unable to describe the elongated structures seen in Fig. |21 

To check that the results are not specific to two dimensions, we have solved the GP 
equation in 3D also. Both ground state configurations are obtained in these simulations. 
As an example we present a density plot of a c. m. state in Fig. |3] Unfortunately, the 
computations are too time consuming for a quantitative mapping of the phase-space. 



V. CORRELATION ANALYSIS 

It is instructive to visualize the dynamics of the many-body wave function Eq. ^ found 
by Wilkin et al. , to see the connection between the wave function and the result of the 
numerical Gross-Pitaevskii calculations (Fig. |21). We therefore calculate the outcome of two 
consecutive measurements, given by the correlation function 

c (.1, z2; tl, t2) = (^U^i,tr)^t(,„t.)^(,„t.)V;(.r,t,)) 

that is, the conditional probability of detecting an atom at (2:2,^2) provided that one was 
detected at (21,^1). Here, we scale the lengths and the energies as we did in Sec. |TT| The 
time is scaled by 1/uj. By making a multi-nomial expansion, the many-body wave function 
(jni) can be expressed in the terms of single-particle harmonic oscillator eigenstates: 

N 

|z,|V2 1 (23) 



^ = E \ Mr, , ' T-^r ■■■^T exp V 
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where the sum is taken over all possible combinations of {rrii} that fulfill mi + ■ ■ - mN 
L. This state can now be easily constructed in the second quantized form |\E') 
f (no, . . . ,11^) \no, . . . ,ni,) where the expansion coefficients are 



1 i_ 

\ 2 



1 / A^l \ 2 / L' 
^^""'•••'"^^ = A^iw^j i (0!r---(L!r J ^'^^ 

and rim denotes the population of atoms in the harmonic oscillator eigenstate of angular 
momentum m. 

When the ideal c. m. state is written as a superposition of different distributions of atoms 
in harmonic oscillator eigenstates, it is natural to determine the correlation by writing 
the field operators in the harmonic oscillator basis 

L ^ 

iP{z,t) = J2 (vrm!)'2 2'"exp(^-|;z|V2)exp(-i(m + l)t)a„. (25) 

m=0 

Here dm is the bosonic particle destruction operator for the state m. The form of the time 
dependent exponential comes from the energy spectrum 

E = uj = \m\+2nr + 1. (26) 

In this case, to minimize the energy, the expansion (f^Hj) contains only terms where the radial 
quantum number Ur is zero. We are interested in the dynamics as a function of the polar 
angle 6 with a fixed distance R from the center of the trap. Hence, we denote zi = Re^^^ 
and Z2 = -Re*^^. In addition we choose ti = 0, 9i = 0, t2 = t and 02 = 0. After doing some 
algebra, one gets 

{i^{z,M)i'{ziM))=T.jl e-^' f {nmiy R'^um (27) 

/ m=0 

for the denominator. The numerator is a bit more tedious but straightforward: 

{lp\zi, ti)lp^{z2, t2)lp{z2, t2)ij{zi, h)) = 

EE E 



L m-1 r) ^2 

\R^'"e-'^' COS ((m - m') (^ - t)) + 



Y^o^,7r^ [m'\{2m-m>)\ (m'ly 

E E E ^^^"^'""^^e-^^ cos iim' -m)ie- t)) + (28) 



f ^ n ,^ m'\m\ 
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2/ IT-m'^m j^2{m'+m) ^-2R'^ , 

f m=Om'>m 



/^^m (^m 1) j^Am^-2B? 

In the first term m" must satisfy certain conditions, namely m" ^ m, m" ^ m' and m" ^ 
m + m' — m" . By combining ()27p and (PHj) . we get a sliglitly complicated expression for 
the correlation function ()22j). To visualize, we plot an example in Fig. El for R = 1. The 
curves illustrate the conditional probabilities at two different values of t when = 6 and 
L = 6. As expected it describes a shape retaining a peaked structure moving clockwise at 
an angular frequency Vt = uj {i. e. f2 = 1 in dimensionless units). 

VI. CONCLUSIONS 

In conclusion, we have studied the ground state properties of an attractively interacting 
zero temperature Bose gas in the presence of an anharmonic trap and a rotational drive. 
The two different ground state configurations, vortex and center of mass rotation were in- 
vestigated at the limit when both anharmonicity and interactions are weak. In this limit, 
analytical conditions for these states to appear were obtained as functions of rotation fre- 
quency, anharmonicity and interaction strength. To go beyond the perturbation theory, 
we have demonstrated that the ground state can be solved by using the ordinary Gross- 
Pitaevskii equation co-rotating with the external drive. This treatment differs from our 



previous work where the coupled equations of motion were obtained by assuming the 
many body wave function to be a product of a cm. wave function and a Bose condensed 
wave function of the motion relative to the cm. Here, we treat the motion of the cm. 
classically. The Gross-Pitaevskii approach enables us to study any strength of the coupling 
or anharmonicity. We note that the comparison between different approximative methods 
shows that the perturbation result is valid for only a very limited parameter range. The 
extension of the perturbative analysis performed in Ref. may improve the situation 
slightly. Finally, we have visualized the cm. dynamics by studying the correlations of two 
consecutive measurements. 

The rotation of a Bose condensate with repulsive interactions in an anharmonic trap is 
already experimentally performed 0|. The condensate was stirred with a laser in a trap 
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FIG. 5: Probability (unnormalized) of observing an atom at 6 when another atom was detected 
at 9 = for t = 0. The sohd hne is the time correlation Eq. (22) at t = and the dashed line 
represents t = 1.0. The unit of time is 1/uj. 

which was created by combining a harmonic magnetic potential with an optical potential 
of Gaussian shape. In principle, by using a similar technique for attractively interacting 



Bose gas one could be able to observe the rotational ground states we have studied. 



possibility to tune scattering lengths with magnetic fields via Feshbach resonances 
offers another tool for probing the parameter space. 
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Since the first submission for publication of the present manuscript, the analysis in the present 
chapter has been extended by Kavoulakis et al. 

Reference jl3l | discusses other types of angular-momentum carrying states; these are, however, 
easily seen to always have higher energy than ^ and ©• 



17 



This figure "figure2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0403326v2 



This figure "figure4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0403326v2 



